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Abstract 

By a suitable transformation, we present the (n + l)-dimensional charged rotating solutions of 
Gauss-Bonnet gravity with a complete set of allowed rotation parameters which are real in the whole 
spacetime. We show that these charged rotating solutions present black hole solutions with two 
inner and outer event horizons, extreme black holes or naked singularities provided the parameters 
of the solutions are chosen suitable. Using the surface terms that make the action well-defined for 
Gauss-Bonnet gravity and the counterterm method for eliminating the divergences in action, we 
compute finite action of the solutions. We compute the conserved and thermodynamical quantities 
through the use of free energy and the counterterm method, and find that the two methods give the 
same results. We also find that these quantities satisfy the first law of thermodynamics. Finally, we 
perform a stability analysis by computing the heat capacity and the determinant of Hessian matrix 
of mass with respect to its thermodynamic variables in both the canonical and the grand-canonical 
ensembles, and show that the system is thermally stable. This is commensurate with the fact that 
there is no Hawking-Page phase transition for black objects with zero curvature horizon. 
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I. INTRODUCTION 



The thermodynamics of asymptotically anti-de Sitter (AdS) black holes, which produces 
an aggregate of ideas from thermodynamics, quantum field theory and classical gravity, 
continues to attract a great deal of attention. One reason for this is the role of AdS/CFT 

n 

duality and in particular, Witten interpretation of the Hawking-Page phase transition 
between thermal AdS and asymptotically AdS black hole as the confinement-deconfinment 
phases of the Yang-Mills theory defined on the asymptotic boundaries of the AdS geometry 
3]. The fact that thermodynamics of black holes may be modified by using other theo- 
ries of gravity with higher derivative (HD) curvature terms provides a strong motivation 
for considering thermodynamics of black holes in HD gravity jjj. The appearance of HD 
gravitational terms can be seen, for example, in the renormalization of quantum field theory 
in curved spacetime |4] , or in the construction of low-energy effective action of string theory 
P). Among the theories of gravity with higher derivative curvature terms, the Gauss-Bonnet 
gravity has some special features in some sense. Indeed, in order to have a ghost-free action, 
the quadratic curvature corrections to the Einstein-Hilbert action should be proportional 
to the Gauss-Bonnet term 0, Q]. This combination also appear naturally in the next-to- 
leading order term of the heterotic string effective action, and plays a fundamental role in 
Chern-Simons gravitational theories ^. From a geometric point of view, the combination 
of the Einstein-Gauss-Bonnet terms constitutes, for five dimensional spacetimes, the most 
general Lagrangian producing second order field equations, as in the four- dimensional grav- 
ity which the Einstein-Hilbert action is the most general Lagrangian producing second order 
field equations 

To investigate the thermodynamics of black holes, one should compute the finite con- 
served and thermodynamic quantities of the system. In general the gravitational action is 
divergent when evaluated on solutions, as is the Hamiltonian and other associated conserved 
quantities. A common approach to evaluating thermodynamic quantities has been to carry 
out all computations relative to some other spacetime that is regarded as the ground state 
for the class of spacetimes of interest. This could be done by taking the original action for 
gravity coupled to matter fields and subtracting from it a reference action, which is a func- 
tional of the induced metric 7 on the boundary dM ^(|. Conserved and/or thermodynamic 
quantities are then computed relative to this boundary, which can then be taken to infinity if 
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desired. Unfortunate 
is not always unique 



y, it suffers from several drawbacks. The choice of reference spacetime 



ll| , nor is it always possible to embed a boundary with a given induced 



metric into the reference background. An extension of this approach which addresses these 
difficulties was developed based on the conjectured AdS/CFT correspondence for asymptotic 
AdS spacetimes |l2,[3|. Applications of AdS/CFT correspondence also include computa- 
tions of conserved quantities for black holes with rotation, NUT charge, various topologies, 
rota ting black strings with zero curvature horizons and rotating higher genus black branes 
[bll Q, Q]- Although the counterterm method applies for the case of a specially infinite 
boundary, it was also employed for the computation of the conserved and thermodynamic 
quantities in the case of a finite boundary [17]. It is believed that appending a counterterm, 
7 ct , to the action which depends only on the intrinsic geometry of the boundary (ies) can 
remove the divergences. This requirement, along with general covariance, implies that these 
terms are functionals of curvature invariants of the induced metric and have no dependence 

n 

on the extrinsic curvature of the boundary (ies) 13]. An algorithmic procedure exists for con- 
structing L± for asymptotic AdS in Einstein gravity and so its determination in this theory 
is unique . But in HD gravity these conterterms have been not introduced till now. Since 
these counterterms should be functions of intrinsic curvature of boundary, the construction 
of these terms for spacetimes with flat boundaries is straightforward \19L l20j[. 

The outline of our paper is as follows. We give a brief review of the field equations of 
Gauss-Bonnet gravity and the counterterm method for calculating conserved quantities in 
Sec. ini In Sec. IIHI we obtain mass, angular momentum, entropy, temperature, charge, 
and electric potential of the (n + l)-dimensional black hole solutions with a complete set of 
rotational parameters and show that these quantities satisfy the first law of thermodynamics. 
We also perform a local stability analysis of the black holes in the canonical and grand 
canonical ensembles. We finish our paper with some concluding remarks. 
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II. THE ACTION AND CONSERVED QUANTITIES 



The gravitational action of an (n+l)-dimensional asymptotically anti-de Sitter spacetimes 
M. with the Gauss-Bonnet term in the present of an electromagnetic field is: 

Ig = -^t f d n+1 x^g~{R -2A + a {R, utJK R^ K - m^RT + R 2 ) - F^F^} 
[ <TaV=7 {e + 2a ( J - 2Q ab e ab ) } 

JdM 

where F^ v = d fl A u — d u A fl is electromagnetic tensor field, A p is the vector potential, A = 
—n(n — l)/2/ 2 is the cosmological constant, a is Gauss-Bonnet coefficient with dimensions 
(length) 2 which we assume to be positive as in heterotic string theory, R puaK) R^ v and R are 
Riemann tensor, Ricci tensor and Ricci scalar of the manifold M. respectively, is induced 
metric on the boundary dA4, 6 is trace of extrinsic curvature of the boundary, Q a bi.l) is 
Einstein tensor calculated on the boundary, and J is trace of: 

J ab = \{Q cd Q cd Qab + 2QQ ac e c b - 2e ac e cd e db - e 2 e afe ) (2) 

The second inte gral is a surface term which is chosen such that the variational principle is 
well-defined |2l|, |22|. In order to obtain the field equations by the variation of the volume 
integral with respect to the fields, one should impose the boundary condition 5A P = on 
dM.. Thus the action (JTJ) is appropriate to study the grand-canonical ensemble with fixed 
electric potential j^J. To study the canonical ensemble with fixed electric charge one should 



impose the boundary condition 5(n a F ab ) = 0, and therefore the gravitational action is 24]: 



Ig = Ig-^- I d n x^n a F ab A b , (3) 



47T 



dMc 



where n a is the normal to the boundary dM. Varying the action or Q with respect 
to the metric tensor g pv and electromagnetic tensor field F^, with appropriate boundary 
condition, the equations of gravitation and electromagnetic fields are obtained as 

Gf, u + Ag^ + aG {2 } = T pv (4) 

V V F^ = (5) 

where G^ v is the Einstein tensor, TjfJ 11 ^ = 2F p ^F pu — \F pu F pu g^ is the energy-momentum 

(2) 

tensor of electromagnetic field and G/J is second order Lovelock tensor: 

G^^ ^^R^uKrRy ^R/ipuaR^ 2R^ ia R v -\- RR^y) 

-\ {R^ K R^ K - AR^RT + R 2 ) g pv (6) 



In general the action Iq, is divergent when evaluated on the solutions, as is the Hamiltonian 
and other associated conserved quantities. Rather than eliminating these divergences by 
incorporating reference term, a counterterm I ct is added to the action which is functional 
only of the boundary curvature invariants. For asymptotically AdS spacetimes with flat 
boundary, all the counterterm containing the curvature invariants of the boundary are zero 
and therefore the counterterm reduces to a volume term as 



ct ~ SttL 



d n xy/^ (7) 
dM 



where L is a length factor that depends on the fundamental constants I and a, that must 
reduce to I as a goes to zero. Thus, the total action, J, can be written as 

I = Ig + la (8) 

n 

Having the total finite action, one can use Brown and York definition 10] to construct a 
divergence free stress-energy tensor as 

T ab = i | (Q ab - 7 ab ) + 2a (3J ab - Jj ab ) + ^-^7 afc j (9) 

To compute the conserved charges of the spacetime, one should choose a spacelike hyper- 
surface B in dM. with metric a%j, and write the boundary metric in ADM form: 

-f ab dx a dx a = -N 2 dt 2 + try (dft + V*dt) (dft + V j dt) , 

where the coordinates <p l are the angular variables parameterizing the hypersurface of con- 
stant r around the origin. When there is a Killing vector field £ on the boundary, then the 
conserved quantities associated with the stress tensors of Eq. (jHJ) can be written as 

Q(0 = / d n - l <P^T ab n a e, (10) 

where a is the determinant of the metric o^- and n a is the unit normal vector to the boundary 
B . For boundaries with timelike (£ = d/dt) and rotational Killing vector fields (£ = d/d<fi), 
we obtain 



M = / eT-VV^n^ 6 , (11) 

JBoo 

J = f d n - l <PV^T ab n% b , (12) 



provided the surface B contains the orbits of (. These quantities are, respectively, the 
conserved mass and angular momentum of the system enclosed by the boundary. Note 
that they will both be dependent on the location of the boundary B in the spacetime, 
although each is independent of the particular choice of foliation B within the boundary 
dAi. In the context of the (A)dS/CFT correspondence, the limit in which the boundary 
B becomes infinite (Boo) is taken, and the counterterm prescription ensures that the action 
and conserved charges are finite. No embedding of the surface B into a reference spacetime 
is required and the quantities which are computed are intrinsic to the spacetimes. 



III. THERMODYNAMICS OF BLACK BRANES 



The rotation group in n + 1 dimensions is SO(n) and therefore the number of independent 
rotation parameters for a localized object is equal to the number of Casimir operators, which 
is [n/2] = k, where [n/2] is the integer part of n/2. Asymptotically AdS solution of the 
field equations (jlj) with cylindrical symmetry with k rotation parameter Oj, can be written 
as 3|: 

/ k \ 2 2 k 

ds 2 = -f(p) I Edt - ^didfa J + ( aidt ~ -fd'Pi) 

\ i=l / i=l 

p 2 k ^ 2 

P ^ f(p) 



k / n 2 

= 2 = E(i + | 



i=l 



. [n - I q . 

M = V 2(n-2)p"- 2 ( » ~ ai ^' (n ° SUm ° n l) ' 

where dX 2 is the Euclidean metric on the (n — k — l)-dimensional submanifold with volume 
LUn-k-i and f(r) is 



™ = 2(n-2Kn-3)a { 1 " f ~ 4( " " 2)( " " 3) ° (? " ^ + ^) >• (13) 

As we will see in the next section, m and q are related to the total mass and total charge 
of spacetimes. Although f(p) for the uncharged solution (q = 0) is real in the whole range 
< p < oo provided a < l 2 /A(n — 2)(n — 3), for charged solution it is real only in the range 
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r o < P < °o where tq is the largest real root of the following equation: 

U4n-3) -£W 2 "~ 2 ' +amr ° ( ""'-"^° < M > 
One may note that the only solution of Eq. ()14j) for the case of uncharged solution is ro = 0. 
In order to restrict the spacetime to the region p > r , we introduce a new radial coordinate 



r as: 

2 2 2 ^ r 2 r 7 

r = p — r =^ dp = — Tjdr 

r z _|_ r ^ 



With this new coordinate, the above metric becomes: 

ds 2 = -f(r) sdt - + [r j 4 To) J2 - zi 2 d&y 



^2 j ^.2 ^,2 ^^,2 

— ^5J(a^0, - a^) 2 + + (r 2 + r 2 )^ 2 , (15) 



, , (r 2 + rg)/(r) 

where now the gauge potential and the metric function are: 



Am = ^ 2(n-2) ( r 2 + r 2 )f- l ((S(5 ° ~ (n ° 011 %) ' (16) 



f( r ) = 777 7W7 on S 1 - \ 1 - A ( n - 2 )( n - 3 ) 7? " , o . ^ + 



2(n-2)(n-3)a P ty^ ' v " ~' v "' (r 2 + r 2 )t ' (r 2 + r 2 )"" lV 

It is notable to mention that the in the static case a, = 0, in contrast with the case of 
uncharged solution, g tt is not equal to g~}. 



A. Properties of the solutions 

One can show that the Kretschmann scalar R^ v xkR^ vXk diverges at r = 0, and therefore 
there is a curvature singularity located at r = 0. Seeking possible black hole solutions, we 
turn to looking for the existence of horizons. As in the case of rotating black hole solutions 
of Einstein gravity, the above metric given by Eqs. ()15 )) -([17 )) has both Killing and event 
horizons. The Killing horizon is a null surface whose null generators are tangent to a Killing 
field. The proof that a stationary black hole event horizon must be a Killing horizon in 
the four-dimensional Einstein gravity [3| cannot obviously be generalized to higher order 
gravity. However the result is true for all known static solutions. Although our solution is 
not static, the Killing vector 

k 

X = d t + YPid^ (18) 
i=i 

7 



is the null generator of the event horizon, where k denotes the number of rotation parameters. 
The event horizon define by the solution of g rr = f(r) = 0. The metric of Eqs. (|15 p -()17 )) 
has two inner and outer event horizons located at r_ and r + , provided the mass parameter 
m is greater than m ext given as: 



n - 1 \ / Tl \ 2 ("-!) n-2 n , . 

m -'- 2 {—2){—2) < 19 ) 
In the case that m = m ext , we will have an extreme black brane. One can obtain the angular 
velocity of the event horizon by analytic continuation of the metric. Setting a« — > iai yields 
the Euclidean section of (fTKj) . whose regularity at r = r + requires that we should identify 
(pi ~ (pi + P+^i, where fij's are the angular velocities of the outer event horizon. One obtains: 

a = ^fr (2°) 
The temperature may be obtained through the use of definition of surface gravity, 



T+ = 27v4 (V ^ )(VMr) ' (21) 



where \ is the Killing vector (jl8|) . One obtains: 



One may note that the temperature (J22J) reduces to the temperature of the uncharged 
solution for the uncharged solutions where r = j^. Next, we calculate the electric 
charge of the solutions. To determine the electric field we should consider the projections of 
the electromagnetic field tensor on special hypersurfaces. The normal to such hypersurfaces 
is 

= i u' = 0, u* = ^, (23) 

and the electric field is = g^ p F pv u v , where N and V 1 are the lapse and shift functions. 
Then the electric charge per unit volume V n -\ can be found by calculating the flux of the 
electric field at infinity, yielding 



The electric potential $, measured at infinity with respect to the horizon, is defined by 
22|: 

$ = I™ - A, X \ =r+ , (25) 



where x is the null generator of the horizon given by Eq. (|18|). We find 



(n-1) 



$ W 2(n-2) s(r 2 +ro2) <^- (26) 
B. Action and conserved quantities 

The action and conserved quantities associated with the spacetime described by (|T3jl can 
be obtained via the counterterm method. Using Eqs. ((TJ), (jSJ) and (J7J), the Euclidean actions 
will be finite provided 

1 _ i Z 2 +4(n-2)(n-3) a-l^/l 2 -4(n-2)(n-3) a , > 

I ~ 3V2(n-2)( n -3) Tj/^T^ 

where we note that L reduces to I as a goes to zero. Denoting the volume of the hypersurface 
boundary at constant t and r by V n -\ = (27r) a; n _fc_i, the Euclidean actions per unit volume 
V n -i in canonical and grand-canonical ensembles can then be obtained through the use of 
Eqs. (0), © and 0. We find: 

P + (r 2 + rg)"- 1 + g 2 P 

167T /2( r 2 +r 2)f-i ' ^ 
J P + (rj. + rp 2 )- 1 - (2n - 3)g 2 Z 2 

167T /2 (r 2 +r 2 )f -l ^ 

We first obtain the action in grand canonical ensemble as a function of the intensive 
quantities (3, ft and $ by using the expression for the temperature, the angular velocity and 
the potential given in Eqs. (J22J), (JUD and (J2HJ): 

, W + ^ ({rl jHL ), (30) 



167r/ 2 (r 2 +r 2 )f~ 1 V (n - 1)(1 - Q 2 Z 2 

where r + in terms of /?, O and $ is 



7r(n - 1) + yVfo - l) 2 - A$ 2 /3 2 (n - 2) 2 ] 
r+ ~ /5AV1 - ft 2 / 2 ' 

Since the Euclidean action is related to the free energy in the grand canonical ensemble, 
the electric charge Q, the angular momenta Jj, the entropy S and the mass M per unit 
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volume V n -i can be found using the familiar thermodynamics relations: 



d$-^y — 2 — g ' (32) 

S = (^- 1 ) I =^ r l + r o)^, (33) 
dl 1 

Ji = -P^t^t = —^nEmai, (34) 

Note that the charged obtained by use of free energy is the same as that obtained in Eq. 
(J2IJ). Of course, if one computes the mass and angular momentum per unit volume V n _i 
through the use of counterterm method introduced in Ref. (explained in Sec. |TlJ), the 
same results will be obtained. 

Equation (|33|) shows that the entropy of the black brane satisfies the so-called area law 
of entropy which states that the black hole entropy equals to one-quarter of horizon area 
. Although the area law of entropy applies to almost all kind of black ho es and black 
strings of Einstein gravity j^, it is not satisfied in higher derivative gravity 3(|. But, for 
our solutions in Gauss-Bonnet gravity with flat horizon the area law is satisfied. 



C. Stability of black branes 

We first obtain the mass M as a function of S, Q and J. Using the expressions (J33J), (J55J) 
and f° r the entropy, the mass and the angular momentum, and the fact that f(r + ) = 0, 
one obtains by simple algebraic manipulation the Smarr-type formula as: 

M(s > J > Q) = ijmky (36) 

where J 2 = | J| 2 = Jf and Z = S 2 is the real positive solution of the following equation: 

(Z - 1)^ - f Mn-l)(n-2)lSJ ] ^ = 

V } 16S 2 \ n[(n - l)(n - 2)S 2 + 2tt 2 Q 2 1 2 ] J ' y } 

It is worthwhile to mention that the thermodynamic quantities calculated above satisfy the 
first law of thermodynamics, 

dM = TdS + n i dJ i + ® d Q- (38) 
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The stability of a thermodynamic system with respect to the small variations of the ther- 
modynamic coordinates, is usually performed by analyzing the behavior of the entropy 
S(M, Q, J) around the equilibrium. The local stability in any ensemble requires that 
S(M, Q, J) be a convex function of their extensive variables or its Legendre transforma- 
tion must be a concave function of their intensive variables. Thus, the local stability can 
in principle be carried out by finding the determinant of the Hessian matrix of S with re- 
spect to its extensive variables, H x . x . = [d 2 S/dXidXj], where X^s are the thermodynamic 
variables of the system. Indeed, the system is locally stable if the determinant of Hessian 
matrix satisfies H^. >x . < 0,13- Also, one can perform the stability analysis through the 
use of the determinant of Hessian matrix of the energy with respect to its thermod yna mic 
variables, and the stability requirement H^. x . < may be rephrased as Hy? y . > [3l|. 

The number of the thermodynamic variables depends on the ensemble which is used. 
In the canonical ensemble, the charge and angular momenta are fixed parameters, and 
therefore the positivity of the heat capacity Cq j = T(dS/dT)q j is sufficient to assure the 
local stability. One can show that the latter criteria is: 

Gr,Q = | (rl + rl)^[ (n (r 2 + r 2 ) M - (n - 2) g 2 Z 2 ) ((r 2 + r 2 )<»"*> + g 2 / 2 ) 
((n - 2)H 2 + 1) \ | (n - 2)q A I 4 ((3n - 6)H 2 - (n - 3)) 

- 2g 2 / 2 (r 2 + r 2 ((3n - 6)S 2 - n 2 + 3) + n(r 2 + r 2 )^ ((n + 2)S 2 - (n + 1)) 

(39) 

Figure Q shows the behavior of the heat capacity as a function of the charge parameter. 
It shows that Cj,q is positive in various dimensions and goes to zero as q approaches its 
critical value (extreme black brane). Thus, the (n + l)-dimensional asymptotically AdS 
charged rotating black brane is locally stable in the canonical ensemble. One may also plot 
T versus S for fixed J and Q, and look for inflection point (s). Indeed, for systems with phase 
transition the dia gram has an inflection point, that demonstrates the region of coexistence 
of various phases [2J|. Since there is no an inflection point in Fig. the black brane is 
stable. 

In the grand-canonical ensemble, the stability analysis can be carried out by calculating 
the determinant of Hessian matrix of the energy with respect to S, Q and J. It is a matter 
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FIG. 2: T versus S for I = 0.09, q = 0.06, a = 0.01, n = 4 (solid line)), n = 5 (dash-dotted), and 
n = 6 (bold line). 

of calculation to show that the black brane is locally stable in the grand-canonical ensemble, 
since the determinant of Hessian matrix, 

f64;2„2 , 64/ 2 , 2 \n-l\ 

Hf Q3 = ^ + ^ (r ° + r+) ] ^- (40) 

bQJ Z 2 ^ 6 (( n _ 2)H 2 + l)(/ 2 g 2 + (r 2 + r 2 )™" 1 )(r 2 + r\ )t n ~ 2 

is positive. It is worthwhile to note that a has not appeared in the thermodynamic quantities 
computed in this section. Thus, the charged solutions in the Gauss-Bonnet gravity has the 
same thermodynamic features as the solutions in the Einstein gravity. This fact is also true 
for the case of static S ph erical l y symmetric blac, hoi. with zer o cu ™t ure _ Q. 
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IV. CONCLUDING REMARKS 



The metric function of charged solution of Gauss-Bonnet gravity introduced in literature 
is not real for the whole spacetime even for static case 3^|. In this paper we presented the 
charged rotating solutions of Gauss-Bonnet gravity which are real in the whole spacetime 
and investigated the thermodynamics of them. We found that in contrast with the case 
of uncharged solution, the temperature is not equal to /'(r + )/47rH. The charged solutions 
may be interpreted as black brane solutions with two inner and outer event horizons for 
m > m ext , extreme black holes for m = m ext or naked singularity for m < m ext , where m ext 
is given in Eq. (|19|) . We found that the Killing vectors are the null generators of the event 
horizon, and therefore the event horizon is a Killing horizon for the stationary solution of the 
Gauss-Bonnet gravity explored in this paper. We computed the finite action of the charged 
rotating solutions through the use of counterterm method. By using the relation between 
the action and free energy, we compute the conserved and thermodynamic quantities of 
the solutions and showed that they satisfy the first law of thermodynamics. Finally, we 
obtained a Smarr-type formula for the mass of the black brane solution as a function of 
the entropy, the charge and the angular momenta of the black brane and studied the phase 
behavior of the (n+ l)-dimensional charged rotating black branes. We showed that there is 
no Hawking-Page phase transition in spite of the charge and angular momenta of the branes. 
Indeed, we calculated the heat capacity and the determinant of the Hessian matrix of the 
mass with respect to S, 3 and Q of the black branes and found that they are positive for 
all the phase space, which means that the brane is stable for all the allowed values of the 
metric parameters discussed in Sec. 11111 It is worth to mention that the charged rotating 
solutions of first order Lovelock gravr^vQ, of second order Loveiock gravity introduced in 
this paper, of third order Lovelock gravity presented in Ref. 20] and the static solutions of 
continued Lovelock gravity |l9| are stable. Thus, we conjecture that the Lovelock terms do 
not have any effects on the stability of black branes with flat horizon. This is commensurate 
with the fact that there is no Hawking-Page transition for a black object whose horizon is 

Q 

diffeomorphic to 1R P and therefore the system is always in the high temperature phase [2]. Of 
course this kind of solutions in the presence of a dilaton field is not stable even in Einstein 
gravity 3- 

Although the thermodynamics of charged rotating solutions of Gauss-Bonnet gravity with 
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flat horizon have been investigated in this paper, no such a kind of solutions with spherical 
horizon have been presented till now. Thus, it is worth constructing these kinds of solutions 
and investigating their thermodynamics. The charged rotating black holes with spherical 

I I 

horizon in Einstein gravity are not stable in the whole phase space [23], and therefore it is 
interesting to investigate the effects of Gauss-Bonnet term on the thermodynamics of such 
kinds of solutions. 
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